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Abstract 

Eta-mesic nucleus or the quasibound nuclear state of an eta {g) meson in a nucleus is caused 
by strong-interaction force alone. This new type of nuclear species, which extends the landscape 
of nuclear physics, has been extensively studied since its prediction in 1986. In this paper, we 
review and analyze in great detail the models of the fundamental 77 -nucleon interaction leading to 
the formation of an ry-mesic nucleus, the methods used in calculating the properties of a bound g, 
and the approaches employed in the interpretation of the pertinent experimental data. In view of 
the successful observation of the 77 -mesic nucleus and other promising experimental results, 

future direction in searching for more 77 -mesic nuclei is suggested. 
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I. INTRODUCTION 


It is well-known that mesons play an important role in nuclear physics. The interaction of 
mesons with nuclei has two complementary components: meson-induced nuclear reactions 
and meson-nucleus bound systems. Thus, an understanding of the ensemble of meson- 
nucleus interactions can enhance our knowledge of nuclear force and nuclear structure. 

The modern era of meson-nucleus physics began with the advent of various meson factories 
in the 1960s, where high-intensity pion (vr) and kaon (K) beams were made available. Since 
then, meson-nucleus bound systems such as vr-mesic and i^-mesic atoms have been studied 
extensively. Consequently, a wealth of information has been obtained about vr-nucleus and 
/f-nucleus interactions. 

For a long time, the role of eta {rj) meson in nuclear physics research was considered 
secondary because the r^-nucleon-nucleon (riNN) coupling constant is much smaller than 
the ttNN and pNN coupling constants. In the mid-1980s, experiments at LAMPF showed 
that 7] mesons are copiously produced in pion-induced nuclear reactions. This led to the 


Q. 


development of the rjN interaction model by Bhalerao and Liu (BL) llH. The model has 
made evident that rj production off a nucleon is dominated by the iV*(1535) resonance and 
that both the rjNN* and ttNN* coupling constants are by no means small. The BL model 
was later used by Haider and Liu to predict the existence of nuclear bound states of the ri 

n 

meson - the ry-mesic nuclei [2[ . 

Mesic nuclei differ from mesic atoms in two imporant aspects. While the formation of 
mesic atoms is driven by the Coulomb interaction between a nucleus and the bound meson, 
the binding of an rj meson into a nuclear orbit is solely due to strong interaction because 
the 7] carries no electric charge. Furthermore, while the size of mesic atoms are of atomic 
scale, the size of ? 7 -nucleus bound systems are of nuclear scale. The prediction of r^-mesic 


nucleus, a novel form of nuclear species adds, therefore, a new dimension to the study of t 


tie 


dynamics of r^-nucleus interaction and the properties of rj meson in nuclear medium |3|, l4|. 

In this paper, we review the progress made in the search of r^-mesic nuclei since its 
prediction in 1986. Among others, we examine the various experimental approaches used in 
the search. We also analyze in depth different methods employed in interpreting the data. 
In section ini we give a comprehensive analysis of the low-energy rjN interaction models that 
are the basis of the formation of r^-mesic nucleus. In section lllli theoretical calculations 
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that led to the prediction of the existence of r^-niesic nuclei are reviewed. In particular, 
we demonstrate the importance of treating realistically the subthreshold r]N interaction 
in a nucleus. Current status on experimental searches for r^-mesic nuclei, including the 
observation of ry-mesic nucleus ^^Mg^, are discussed in section HVl Suggestions on future 
search for ry-mesic nuclei are given in section |Vl 


II. LOW-ENERGY ETA-NUCLEON INTERACTION 


The threshold of r^-nucleon system is 1488 MeV which is 47 MeV below the S'!! baryon 
resonance TV* (1535). This resonance couples strongly to the r]N system with an 7]N decay 
branching fraction of about 45-60%. Consequently, in the threshold region the rjN inter¬ 
action is dominated by iV*(1535) and it is attractive. Bhalerao and Liu [l| formulated an 
off-shell isobar model for threshold pionic rj production on a nucleon and rjN scattering. They 
treated the three dominant reaction channels - ttN, ttttN, and rjN - in a coupled-channels 
formalism and unitarized the model through the generation of the coupled T-matrices. The 
parameters of the model were determined from htting only the ttN phase shifts and inelas¬ 
ticity parameters in the -P 33 , -Pii, and S'!! channels over a broad range of energies. With these 
determined parameters, the model was used to predict the 7r~ +p ^ 7] +n cross sections and 
the rjN scattering length, ar,N- The predicted scattering length has a positive real part.^ 

Because the rjN interaction is attractive, the inequality Re[a^Ar] > 0 indicates that there 
is no s-wave r^-nucleon bound state a. Howeve., it can have an interesting nuclea.* 
tion. As will be shown in the next section, a first-order r^-nucleus optical potential, 14^^, is 
proportional to with trjN being the f-matrix of the rjN scattering and Fa the nuclear 

form factor. Consequently, Re[a^Ar] > 0 leads to Re[f^jv] < 0 and, thus, to Re[f4yi] < 0, 
he., to an attractive ? 7 -nucleus interaction. The attraction, if strong enough, opens the 
possibility of having an rj bound in a nucleus to form a short-lived r^-nucleus bound state. 
Indeed, the first prediction of the existence of such nuclear bound states, the u-mesic nuclei, 
was made by Haider and Liu 

The pN scattering length has since been extensively studied by many researchers. In 


Table [T] we list some representative published results ij, I6l-l27l| . Owing to the unavailibity of 


[1] The sign convention ft{p) +aip'^^, where / denotes the rjN scattering amplitude, was used. 
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TABLE I: Eta-nucleon s-wave scattering lengths ar^N- 


ar^N [fm] 

Method Used“ 

OSE 

Reference 

0.20-bz0.26 

ChPT with pseudo-potential 

-b-b 

Kaiser et al. [6] 

u.zia_g Qgg -|- iu.zoo_Q Q55 

Chiral unitary model 

-b 

Mai et al. [7] 

0.26-bi0.25 

Chiral unitary approach 

-b 

Inoue et al. [8] 

0.27-bi0.22 

Coupled-channel Isobars 

-b-b 

Bhalerao and Liu [Ij 

0.28-bz0.19 

Coupled-channel Isobars 

-b-b 

Bhalerao and Liu [Ij 

0.30-bi0.18 

Coupled-channel T-matrices 

-b-b 

Durand et al. [9j 

0.32-bz0.25 

Chiral EL 

-b 

Ramon et al. [10] 

0.3781°:?^?+i0.20llO;0^3 

ChPT 

-b 

Mai et al. [7] 

0.41 -bzO.26 

MEM 

-b-b 

Gasparyan et al. [llj 

0.41 -biO.56 

A-matrix, solution G380 

-b 

Arndt et al. [12] 

0.42 -b iO.34 

MEM 

-b 

Sibirtsev et al. [13] 

(0.476 - 0.481) -b i(0.279 - 0.289) 

Final-state interaction 

- 

Faldt and Wilkin [1^ 

0.487-biO.171 

EL / A-matrix 

- 

Feuster and Mosel [151 

0.51 -bz0.21 

EL / A-matrix 

-b 

Sauermann et al. [1^ 

0.52-bi0.25 

Final-state interaction 

- 

Willis et al. [17| 

0.54 -b iO.49 

ChPT 

-b 

Krippa [1^ 

(0.55 ± 0.20) -b i0.30 

Final-state interaction 

- 

Wilkin [1^ 

0.577-biO.216 

A-matrix 

- 

Feuster and Mosel [151 

0.68 -b iO.24 

ChPT with pseudo-potential 

-b-b 

Kaiser et al. [201 

(0.710 ± 0.030) -b i(0.263 ± 0.025) 

Coupled-channel T-matrices 

-b 

Batinic al. 

0.75(4) -b i0.27(3) 

A-matrix 

-b 

Green and Wycech [22] 

0.772(5) -bi0.217(3) 

ChPT 

-b 

Nieves and Arriola [231 

0.87-bi0.27 

A-matrix 

-b 

Green and Wycech [24] 

0.91(6) -bi0.27(2) 

A-matrix 

-b 

Green and Wycech [2^ 

0.98-bi0.37 

Quark model 

-b-b 

Arima et al. [2^ 

1.03-bz0.49 

Vector-meson dominance 

-b 

Lutz et al. [2^ 

1.05-bz0.27 

A-matrix 

-b 

Green and Wycech [24] 

1.14-bi0.31 

A-matrix, solution Fit A 

-b 

Arndt et al. [12] 


“ OSE: Off-shell Extension (see the text); ChPT: Chrial Perturbation Theory; EL: Effective Lagrangian; 
MEM: Meson-exchange Model. 


an T] beam, cannot be extracted directly from r^-nucleus to r^-nucleus elastic-scattering 
experiments; rather it has to be inferred from experimental data having an rj in the final 
state by way of using theoretical models. 

An inspection of Table H] shows that the values of Im[a^ 7 v] are confined into a narrower 
range than that of Re[a,y 7 v] which varies between 0.20 and 1.14 fm. The strong model 
dependence of Re[aj,Ar] is due to the fact that Re[a^ 7 v] is not directly constrained by the 
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data. This dependence is further evidenced by the noted large differences between the 


various calculated scattering lengt 
varies from 0.87 to 1.05 fm in Ref. 
and Gridnev 


IS given in some same publications. For example, Re[a^Ar] 


24| . To emphasize the strong model dependence, Birbrair 


28l | calculated using two reaction mechanisms. They showed that one 


of the mechanisms, the iV*(1535) resonance mechanism, gave Re[a,,Ar] = 0.56 fm. The 
other mechanism, the oo-meson-exchange mechanism, gave Re[ar;Ar] = —0.15 fm (even the 
sign changed). The combined effect of the two mechanisms on the real part of the scattering 
length was, however, not given. We note that the oo-meson-exchange has also been included 


in some meson-exchange models (MEM) 


11 


13| , although the individual contribution of the 


exchange diagram was not mentioned. In view of the result of MEM, we believe that an 
overall negative Re[ar,Ar] is unlikely. 

The second column of Table |T] indicates the model or method used in determining the 
scattering length. The possibility of making an off-shell extension of the model/method is 
given in the third column of the table. An off-shell extension (OSE) of the r]N model is 
particularly important for the investigation of ry-nucleus interaction. This is because for the 
formation of an ? 7 -nuclear bound state, the basic 7]N interaction is off-shell. The models 
having both the off-shell momentum and off-shell energy dependences are indicated with a 
double plus sign (-I--I-) in the OSE column. Many models do not have off-shell momentum 
form factors. However, as the lack of an explicit off-shell momentum form factor is equivalent 
to an off-shell momentum form factor having an inhnite range in the momentum space, these 
models are labeled with a single plus sign (-I-) so long as they can be used to calculate rjN 
interaction at subthreshold energies. Otherwise, the models are labeled with a minus sign 
(-), as is the case with the final-state-interaction model Q,Q- 

As can be seen from Table [H many models are based on the iF-matrix approach [l2, 15, 


22 


24 |25| . A discussion on this approach is, therefore, in order. Within the context of 


the JF-matrix approach, one often begins by parametrizing the K matrix and then relates 
it to the T matrix by Heitler’s damping equation T = K — mK5{E — Ho)T, where Hq is 


the free Hamiltonian 


29|. One then uses the T matrix to fit the data, whence to determine 


the parameters of the model. However, the data fitting can determine only the on-shell T 
matrix and hence, by way of Heitler equation, only the on-shell K matrix. In other words, rjN 
models formulated within the framework of iF-matrix approach do not provide an explicit 
way of making off-momentum-shell extension. The Heitler on-shell relation has also led to 
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the caution on the uniqueness of the off-shell T matrix obtained by hrst extrapolating the 
on-shell K matrix (determined from fitting data) to off-shell (subthreshold) energy region 
and then to infer the subthreshold T matrix from the extrapolated K matrix 30j. 


Besides it'-matrix and T-matrix methods, many ant 


ory (ChPT) to calculate the rjN scattering length 


3 , 


lors applied chiral pertubation the- 


18 |, 


20 


23l |. The results are also given 


in Table [H The corresponding Re[a,^Ar] vary between 0.20 to 0.77 fm, depending on how the 
leading orders in the chiral expansion were calculated. However, as pointed out by Kaiser 
et al. 1^, in order to investigate the formation of resonances one needs a non-perturbative 
approach which sums a set of diagrams to all orders. This summation is beyond the frame¬ 
work of systematic expansion scheme of ChPT. To overcome this difficulty, a combination 


of ChPT and pseudopotential methods was employed in Refs. I d l20l . It is interesting to 


note that the model of Ref.j^ is an improvement of that in Ref. 


01 


and it contains more 


terms. This improved model leads to a smaller 7]N scattering length; = 0.20 -|- i0.26fm. 
The authors of Ref. believe that the smaller scattering length is a result of cancellations 
among the various reaction diagrams in the model. 

It is further noteworthy that in a recent coupled-channel isobars approach, Durand et 
al. included five meson-baryon channels and nine isobar resonances. By fitting directly 
the 7r~p —?■ rjn data, their model gave an = 0.30 -|- iO.18 fm, which is remarkably close 


Q 


to that given by the BL [Ij model.^ 

In summary, there are compensations among contributions from various reaction mech¬ 
anisms considered in the riN models. These mutual cancellations could be the reason that 

fi n 

the two most recent models [6|, l9[, which contains many reaction diagrams, give rise to a 
small QriN- In this respect, we surmise that the BL model has grasped the essence of the 
r]N dynamics. Clearly, the quality of the rjN models will be ultimately determined by their 
predictive power of the binding energies and widths of r^-mesic nuclei. 


III. ETA-NUCLEUS INTERACTION AND ETA-MESIC NUCLEUS 

The wave function of an r^-nucleus bound state (or an r^-mesic nuclear state) \l/„ satishes 
the eigenvalue equation where 77 = "Ho + V is the Hamiltonian and Sn is 

[2] The BL model fits the ttN phase shifts and inelasticities instead of the tt~p —>■ rjn data. 
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the eigenenergy. (For simplicity of notation, the quantum label n will be omitted, but 
understood). Solving the above equation is equivalent to solving the integral equation = 
Gq where Gq is the free Green’s function. We will discuss and analyze the physics 
contents of three different approaches used to solve the eigenvalue equation in order to 


calculate the eigenenergy £ 




A. Covariant eta—nucleus optical potential 


Solving a four-dimensional eigenvalue equation requires a full relativistic description of 


the nucleus, which is still not available at this time. Consequent 


reduction 


32l | to obtain a covariant three-dimensional equation 


y, we make a covariant 


3ll: 


— M) + / dk<k' \V\k> ^p(k) = £^p(k') . 

2/i J 


(3.1.1) 


Here, k, k', and /i are, respectively, the initial, hnal relative momentum, and the reduced 
mass of the r^-nucleus system. We denote the eigenenergy £as£ = E — iT/2 with i?(< 0) 
and r(> 0) representing, respectively, the binding energy and width of the //-nucleus bound 
state. In spite of its Schrodinger-like form, Eq. fl3.1.ip is fully covariant. The main advantage 
of working with a covariant theory is that the //-nucleus interaction V can be related to the 


elementary r]N process by unambiguous kinematical transformations 33l. l34 1 


The three-dimensional covariant matrix elements < k' | 1/ | k > in Eq. fl3.1.ip is related 
to the fully relativistic one by 


where 


and 


<k' \V\k>= V^(k'2) < k' I V{W, k^)\k> y/^(k2) , 

VF = -y/Tf 2 -)- ^ , 

Mrj + Ma 


R{k^) = 


E^{k) + EA{k) 


(3.1.2) 


(3.1.3) 


(3.1.4) 


In Eq. fl3.1.3p . Kj. is the magnitude of the on-shell //-nucleus relative momentum. At the 
//-nucleus threshold, Kr=0. 
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The three-dimensional relativistic wave function ip is related to the fully relativistic one 


by 

^(k) = ^(k. • (3-1-5) 

As a result of the application of the covariant reduction, the zeroth components of the 
four-momenta k and k' are no longer independent variables but are constrained by 

k° = W-EA{k), k'^ = W - EA{k') . (3.1.6) 

The hrst-order microscopic r^-nucleus optical potential is represented by the diagram 
shown in Fig{T] It can be expressed in terms of the r]N interaction, namely. 


< k' I V{W) I k > 


= ^ MQ < k', -(k' -h Q) I t{y/^)rjN^r,N \ k, -(k Q) > 

j 

X 0*(-k'-Q)0,(-k-Q) , (3.1.7) 


where the off-shell riN interaction trjN^TjN is weighted by the product of the nuclear wave 
functions <p*j(f>j corresponding to having the nucleon j at the momenta —(k-|-Q) and —(k'-|-Q) 
before and after its collision with the rj meson, respectively. The is the r]N invariant 
mass and is equal to the total energy in the c.m. frame of the rj and the nucleon j. It is 
given by 


s, = [{lT-Ee,,(Q)r-Ql 




Mjj Mj\f — I 6j 


2Mc,j 


/ + Ma 

\Mri T Afjv 


2 

< (M, + , 


(3.1.8) 


where \ej\ is the seperation energy of nucleon j. The Q, and Mc,j are, respectively, 

the momentum, total energy, and mass of the core nucleus arising from removing a nucleon 
j of mometum — (k-1- Q) from the target nucleus of momentum —k. At the threshold of the 
//-nucleus system, W = -\- Ma- 

Equation (13.1.71) indicates that the calculation of V involves integration over the Fermi 
motion variable Q and, hence, the matrix elements of t^N^r^N are to be calculated at both the 
on-shell and off-shell momenta . On the other hand, Eq. (I3.1.8p indicates that the evaluation 
of niust be Carried out at subthreshold energies. 












k, W - EA{k) k', W - EA(k') 



FIG. 1: Diagrammatic representation of the process which determines the first-order covariant rj- 
nucleus optical potential. The dashed line is an rj, the solid circle is the r]N scattering amplitude, 
and the open circles are nuclear vertex functions. The heavy lines represent nuclei which are kept 
on their mass shell (denoted by a cross). 


The off-shell matrix element of tnN^rjN in the r^-nucleus system is related to the off-shell 
riN scattering amplitnde A in the rjN system by 


<k',-(k' + Q) 


\ k, (k+Q)> 


V Er,{p')EN{p')Er,{p)EN{p) 
V^,(k0^iv(k' + Q)E,(k)E^(k + Q) 




(3.1.9) 


where p and p' are the initial and final relative three-momenta in the c.m. frame of the rjN 
system. As already mentioned, the kinematical transformations between the variables on 
the left-side and right-side of Eq. fl3.1.9p are nnambignons in a three-dimensional covariant 
theory. 

We define the on-shell limit as |p'| = |p| = Po and = Erj{po) + Ej^{po) = where 
Po is the on-shell pN relative momentnm. A natnral way of parameterizing A is 


-^(\/^,p',p) = - 


47rV^r,(p')^Ar(p')^r,(p)^iv(p) 




(3.1.10) 


so that in the on-shell limit {da/dVL)r,N^riN =| P- The E has the standard partial-wave 
expansion of a spin 0-spin 1/2 system: 

•T(^,p',p) = + (^+l)4r,2j+(\/^)P')P) ) 

- ia- (p X p') {tiT,2jAV^E',p)-tiT,2j+i\/^^P'^P) ) P'A^) ] ) (3.1.11) 


where p' = |p'|, p = |p|, z = p ■ p', = £ ± 1/2, and T is the isospin of the pN system and 

eqnals to 1/2. In the on-shell limit, 
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4T,2i±(\/^,P',P) 


(exp[ 0 


(3.1.12) 


Vp ^ ' 2*^0 

The phase shifts 6^ are complex-valued because the thresholds for pN ttN and pN httN 

reactions are lower than the threshold for rjN scattering. When p —)■ 0, <5^ —)■ and 

n.,, f, /T7. O. o') 

(3.1.13) 


^2T,2j±i.y/^^P^P) 21 (£) 

)■ p a^Taj± 


P 

The 02^ 2 j a^ 2 T 2 j respectively, the (complex) pN scattering length and volume. Near 

the threshold, only the s-wave term, in Eq. fl3.1.1ip . is important. 

Different off-shell models give different off-shell extensions of A to kinematic regions 
where p ^ p' and ^ y/Wo . In the seperable model of Ref.[l|, the off-shell amplitude is 
given by 


with 


ta{y/^,p',p) - K{y/^,P',P) ) ’ 

(3.1.14) 

K - jEr,{p')EN (pO Erj (p) En (p) , 

(3.1.15) 

Na = K{^j,p')K{^j,p) OC {p'pYv^{p')Ve{p) , 

(3.1.16) 

= yi- - - E;(y3-) - e;( V37) - E;,(y3-), 

(3.1.17) 


and 


Here a is a short-hand notation for the quantum numbers (£, 2T, 2j) of the isobar resonance. 
The Mq, is the bare mass of the isobar a and E", E", and E“^ in Eq. fl3.1.17p are the self¬ 
energies of the isobar associated, respectively, with its coupling to the pN, ttN, and titiN 
channels |l|. The coupling constants and form factors are denoted by g and v. At the pN 
threshold, only the s-wave pN interaction is important, which limits the isobar to A^*(1535) 
or « = (£,2T,2j) = (0,1,1). 


The full off-shell calculation of S was carried out using the inverse-iteration method 351 

n 

and the pN model of Bhalerao an^ Liu [1[. The calculation showed that the existence of 


p-mesic nuclei is indeed possiible j2|. This possibility was reaffirmed by Li et al. 3^ who 


employed a different method, the Green’s function method. The results obtained with the 
inverse-iteration method and with improved numerical integration techniques over the Fermi 
motion variable Q of the nucleon are given in Table UTl No bound state solutions were found 
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TABLE 

off-shell 


II: Binding 
calculation 


energies and half-widths (both in MeV) of ry-mesic nuclei given by the full 


3ll |. No bound-state solutions were found for mass number A < 12. 


Nucleus 

Orbital (ni) 

E-iT/2 

12c 

Is 

-(1.19-hi3.67) 

160 

Is 

-(3.45 -Fi5.38) 

26 Mg 

Is 

-(6.39-Fi6.60) 

40 Ca 

Is 

-(8.91 -Fz6.80) 

90 Zr 

Is 

-(14.80-Fz8.87) 


Ip 

-(4.75 -Fi6.70) 

208 pb 

Is 

-(18.46 -FilO.ll) 


2s 

-(2.37-hi5.82) 


Ip 

-(12.28-hi9.28) 


Id 

-(3.99-hi6.90) 


for nuclei with mass number A < 12. The ^stematic feature of having more bound states 
in heavier nuclei has been discussed in Ref.|2j. In short, it is the increasing compactness of 
the nuclear system with the mass number as well as the increasing magnitude of Re[a,,^ 
that help in the formation of an p-mesic nucleus. In particular, so long as the BL model I| 
is used, no bound state is possible in nuclei lighter than 


B. Factorization of covariant optical potential 


Once the full dynamics of the r^-nucleus optical potential has been understood, it is of 
interest to see whether more insight could be gained from using a simplihed theoretical 
formalism. In this respect, a factorization approximation (FA) has been proposed by Liu 


and Haider 


31| . Within the context of FA, the ijN scattering amplitude in Eq. fl3.1.7p is 


taken out of the Q-integration at an ad-hoc hxed momentum (Q) and an ad-hoc energy \fs\ 
< k' I F(V^) I k >=< k', -(k' + (Q)) I I k, -(k + (Q)) > /(k' - k) , (3.2.1) 

where 

/(k' - k) = - Q),#.,(-k - Q) (3.2.2) 
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is the nuclear form factor having the normalization /(O) = A. In Eq. fl3.2.ip . the off-shell 
is still defined by the same functional dependences on various momenta and energies 
as given by Eqs. fl3.1.9p and 03.1.141) . except that Q and are now replaced by (Q) and 
■y/i, respectively. The choice of (Q) is certainly not unique. It was suggested in Ref. |3l| to 
take an average of two geometries corresponding, respectively, to having a motionless target 
nucleon fixed before and after the rfM collision. This leads one to set 


(Q) = - (^) (k'+ k) , (3.2.3) 

This choice has, in addition, the virtue of preserving the symmetry of the t-matrix with 
respect to the interchange of k and k'. Because Eq. 03.1.8p shows that the rjN interaction in 
a nucleus occurs at subthreshold energies, it is therefore reasonable to set 

Vl = Mr, + Mn - A , (3.2.4) 

with A being a phenomenological energy-shift parameter. From eq. fl3.1.8l) one sees that 


A = 




Mr, -|- Ma 


2Mcj \A^r) A M]\f 


= (Bn) , 


(3.2.5) 


where the average, denoted by ( ), is over all the nucleons (j = 1,...,A). The A or {B^) has the 
meaning of averaged binding of the target nucleons. It is worth noting that the subthreshold 
nature of the hadron-nucleon interaction, Eqs. (l3.1.81 and (I3.2.4p . is also evident in Ji'-mesic 


and TT-mesic atoms. We refer to section III of Ref. 


3l| for details. 


In Table IIIIl we present the bound-state solutions obtained from using the factorized co¬ 
variant potential V [Eq. fl3.2.ip ] with A = 0,10, 20, 30 MeV, and with all the other interaction 
parameters being the same as those for Table [TTl The nuclear form factors used in the calcu¬ 


lations are from Refs. 


31 


U, l37j. A comparison between Tables UTI and IIIIl indicates that the FA 
results obtained with A = 30 MeV are very close to the full dynamical calculation results. 
This value of A is similar to the one found in pion-nucleus elastic scattering studies [38 1. 
From the results in Table m and from Eq. fl3.2.4l) we conclude that the r]N interaction in a 
nucleus takes place mainly at an energy about 30 MeV below the rjN threshold, Mr, -|- Mjv. 
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TABLE III: Binding energies and half-widths of ry-mesic nuclei obtained with the factorization 
approximation for different energy-shift parameters A (in MeV). There is no bound state in nuclei 
lighter than 


Nucleus 

Orbital {ntj 


E-iT/2 

(MeV) 


A = 0 

A = 10 

A = 20 

A = 30 

12 c 

Is 

-(2.18-bi9.96) 

-(1.80-bi6.80) 

-(1.42-bi5.19) 

-(1.10-bi4.10) 

160 

Is 

-(4.61 -bill.57) 

-(3.92 -bz8.13) 

-(3.33-bi6.37) 

-(2.84-bi5.17) 

26Mg 

Is 

-(10.21 -bil5.41) 

-(8.95-bill.17) 

-(7.94-bi8.97) 

-(7.11 -bi7.46) 

"OCa 

Is 

-(14.34 -bil7.06) 

-(12.75 -bil2.55) 

-(11.53-bil0.21) 

-(10.51 -bi8.59) 

90Zr 

Is 

-(21.32 -bil8.59) 

-(19.15 -bil3.97) 

-(17.58-bilL54) 

-(16.29-bi9.84) 


Ip 

-(8.27-bil6.01) 

-(7.19-bill.47) 

-(6.23-bi9.48) 

-(5.40-bi7.94) 

208pb 

Is 

-(24.06 -bil9.18) 

-(21.88-bil4.44) 

-(20.28-bilL96) 

-(18.96-bilO.22) 


2 s 

-(4.89-bill.04) 

-(3.67-bi8.28) 

-(2.81 -bi6.79) 

-(2.12-bi5.72) 


Ip 

-(18.33 -bil8.97) 

-(16.31 -bil4.27) 

-(14.81 -bill.79) 

-(13.56-bil0.06) 


Id 

-(8.27-bil4.07) 

-(6.17-bil0.56) 

-(5.58-bi8.71) 

-(4.66-bi7.41) 


C. Static approximation 


One special case of factorized potential is the static approximation to the potential. In 
the static approximation, not only is the ryiV amplitude factorized out of the integration in 
Eq. fl3.1.7p . but also all the hadron masses are treated as being massive with respect to the 
momenta. The static approximation was first used in the study of mesic atoms 39|, where 
the isospin-averaged spin-nonflip part of the hrst-order static optical potential for a spin-0 
hadron has the general form 40 | 


< k' I I k > = - 


X 


E 

^=0,1 


dvr^/i 
(2 £+l)| k' 


1 + 


rrih 

niN 


(1 -b rrih/mN) 
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f) 

‘‘hN 


/(k' - k) 


P.(k' ■ k) , 


(3.3.1) 


where mh is the hadron mass, /i the hadron-nucleus reduced mass, the effective £-th 
partial-wave hadron-nucleon {hN) amplitude, and Pi is the Legendre polynomial of order 
i. In Eq. fl3.3.ip . k and k' are the initial and hnal hadron-nucleus relative momenta. It is 
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instructive to see the relation between this last equation and the fully covariant amplitude, 
Eqs. fl3.1.9p - fl3.1.1ip . We hrst note that in the static approximation the target nucleon is 
treated as being at rest before as well as after its collision with the hadron 39|. Hence, the 
initial and hnal hN relative momenta, p and p', in the c.m. frame of the hN system are 
p = k/(l + mh/misf) and p' = kY(l + mh/misf), respectively. Clearly, p' • p = k' • k. In 
terms of the variables p and p', Eq. fl3.3.ip becomes 

1 


< k' I I k > = - 


1 + 


dvr^/i 

X 5^(2£ + 1)|p'|'|p 
^= 0,1 


rrih 

rriN 

‘'hN 


/(k' - k) 

■ P) . 


(3.3.2) 


For r^-mesic nuclei calculations, rrih = mri, £ = 0, and Hence, 


It is easy to see that when mass of the z-th particle is treated as being very massive with 
respect to its momentum such that Ej ~ ruj, the multiplicative factor in front of the am¬ 
plitude A in Eq. fl3.1.9p becomes unity and that in front of the amplitude in Eq. fl3.1.10p 
becomes 


rriN + rur, 


An'^mN'mri dTr^m^ 


1 + 


rur. 


rriN 




1 + 




rriN 


(3.3.4) 


which is exactly the multiplicative factor in Eqs. fl3.3.ip - fl3.3.3p . 

In mesic-atom studies, the on-shell hadron-nucleon scattering length was often used 
for o/iTv in Eq. (l3.3.ip . We call the use of on-shell scattering length the on-shell static ap¬ 
proximation. It was pointed out by Kwon and Tabakin |35j that atN should be regarded as 
an effective amplitude. In what follows, we will denote this effective amplitude as while 
use QrjN to denote exclusively the on-shell scattering length. 

We show in Table IIVI the calculated results given by the static apmoximation (SA) for 
light nuclei. The results of using a^Ar= 0.23 -|- i0.09 fm (BL model ij] at 30 MeV below 
threshold), a^Ar= 0.28 -|- i0.19 fm (BL model at threshold), and arjN= 0.48 -|- i0.08 fm 
(from Fig.l of Green-Wychech (GW) j^| model at 30 MeV below threshold) are given, 
respectively, in the third, fourth, and hfth columns. We have found that the use of arjN= 
0.91 -l-i0.27 fm (GW model at threshold) can give bound state in ^He (result not shown). 
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TABLE IV: Binding energies and half-widths of ? 7 -mesic nuclei (Is orbital) obtained with static 
approximation for different values of on-shell scattering length and effective off-shell amplitude 
UriN (in fm). A dash (-) indicates the absence of bound state. 


Nucleus 

Nuclear Form Factor [37] 

E-ir/2 (MeV) 

BL model 

GW model 

o-rjN = 0.23 -|- z0.09 

0.r]N = 0.28 -|- f0.19 

^riN = 0.48 -|- z0.08 

^He 

Hollow exponential 

— 

- 

- 

^He 

3-parameter Fermi 

- 

- 

-(6.02 -Fi3.37) 

^Li 

Modified harmonic well 

- 

— 

-(3.58-Fi2.05) 

^Be 

Harmonic well 

- 

- 

-(12.55-hf3.72) 

10b 

Harmonic well 

-(0.50-hi2.72) 

-(0.93-Hf8.70) 

-(14.37-Hf3.84) 

12c 

Harmonic well 

-(1.71 -Fi3.51) 

-(2.91 -hilO.22) 

-{17.71 +am) 

160 

Harmonic well 

-(3.44-Fz4.24) 

-(5.42-FZ11.43) 

-(21.02-Ff4.19) 

26Mg 

2-parameter Fermi 

-(7.75 -Fz5.89) 

-(11.24-FH4.76) 

-(30.07-Ff4.89) 


However, owing to the rapid decrease of rfM amplitude of the GW model with the energy, 
one sees from Table CYl that there is no bound state in ^He when ar/N of the GW model is 
used. 

For both the BL and GW models, a^jN < CLr^N- In fact, the decrease of the rjN interaction 
strength at subthreshold energies is a very general feature. Furthermore, this decrease is 
very model dependent. Hence, it becomes impractical to employ the models listed in Table |T] 
for which the off-shell dependence cannot be easily reconstructed from the corresponding 
publications. 

Upon comparing columns 3 and 4 of Table HVl with the columns of A=30 and A=0 MeV 
of Table IHIl respectively, we see that the corresponding binding energies and half-widths 
are quite similar to each other. This is to be expected because, as discussed above, SA is 
obtained from the factorization approximation in the infinite-mass limit. Quantitatively, SA 
gives slightly stronger binding energies. This is because in SA there is no off-shell momentum 
form factor at the rjNN* vertices. This slight difference between the SA and FA causes a 
notable difference for the “boderline” nucleus. For example, while the SA predicts a loosely 
bound f] in the FA predicts no bound state in ^°B. Hence, one has to exercise caution 
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TABLE V: Values of for nuclei having a mass number A <10. 


Nucleus 

Nuclear Form Factor [3^ 

a™" (fm) 

^He 

Hollow exponential 

0.49 -h i0.09 

^He 

3-parameter Fermi 

0.35 i0.09 

^Li 

Modified harmonic well 

0.35 -h i0.09 

^Be 

Harmonic well 

0.24 i0.09 


in interpreting the calculated results, particularly if the results indicate a loosely bound 
ry-mesic nucleus. 

It is, thus, informative to determine the smallest value of an effective rjN scattering 
amplitude, o^"', that can bind the 7] into the Is nuclear orbital in a given nucleus. Although 
the real and imaginary parts of this minimal off-shell amplitude are not independent of each 
other, we set, without loss of generality, Im[a,,jv] to be 0.09 fm, as suggested by the two 
off-shell models (BL and GW) discussed above. We then searched for Re\a^^]. The results 
are given in Table |V] for nuclei having mass number A=3 to 9. We choose this mass range 
because of the existing strong interest in finding r^-nuclear bound states in light nuclear 
systems. 

In the literature, the static approximation is sometimes referred to as the local-density 
approximation (LDA) in which the nucleus is treated as an infinite and motionless nuclear 
matter. Since the pioneering study of pionic atoms by Ericson and Ericson [^ . the LDA 
has been extensively applied to studying the vr", K~, and p-atoms. In recent years, the 
method has been applied to the investigation o 


' p-mesic nucleus in and heavier nuclei 
. One calculates the binding energy of the 
? 7 , if it exists, by solving numerically 48l-l45| the coordinate-space Klein-Gordon equation 


by Garcia-Recio et al. 4]| and Gieply et al. 42 


+ Urj{Re[Urj], p{r))] 'll) = , 


(3.3.5) 


where n^(Re[w,j], p(r)) is the self-energy of the rj. In Eq. (l3.3.5p Urj denotes the total energy 
of the 7] and is defined by uJri = {Bn + mn) — iT/2, where Bn is the binding energy of the 
r]. The self-energy 11^ is related to the optical potential R? by n,^(i?e[aJ^],p) = 2i?e[aJ^]14,. 
Theo.etical .nodel. fo. used in Refs.Q and Q aae very diffeept. leading quite 
different results. One notes, among others, the models used in Ref.|4]| gave rise to very 
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large widths while the models in Ref. ^ gave narrow F. A common hnding in both these 
references is the strong sensitivity of the calculated to the energies at which the r]N 
interaction takes place. Their hnding agrees with our discussion of Tables IIIII to |Vl 

In summary, it is important to use an effective off-shell r]N amplitude at the appropriate 
subthreshold energy. Since the energy dependence of the off-shell amplitude is highly model- 
dependent, experimental determination of the lightest nucleus in which an rj can be bound 
constitutes one of the many ways to differentiate various theoretical r]N models. 


IV. EXPERIMENTAL SEARCH FOR 77 -MESIC NUCLEI 


The unavailability of an r^-meson beam makes the hadron-induced nuclear rj production 


the sole way to study r^-nucleus interaction, including the 
Production of rj by pions [4^, protons |47H50| . and deuterons 
in various laboratories. 


brmation of r^-mesic nuclei. 


51 


53l have been carried out 


Experiments in search for ? 7 -mesic nuclei can be divided into two types. In the hrst type of 
experiment one looks for a peak in the spectrum of an emerging particle as a signature of the 
formation of an r^-mesic nucleus. In the second type of experiment the bound state of rj is not 
measured and the experimental hnal state is composed in fact of an unbound rj and a nucleus. 
The total and differential r] production cross sections are measured as a function of the r]- 


5^ to 


nucleus relative momentum One then uses the Watson hnal-state-interaction method 
infer from the data whether there is an ? 7 -nucleus bound state. 

In this section, we will discuss in detail some representative experiments in each of the 
above categories. Our emphasis is on the methods used for theoretical analysis of the data. 
For detailed technical aspects of the experiments, we refer the readers to the excellent review 
by Machner jd]. 


A. Spectral methods 


The hrst search for r^-mesic nucleus was carried out by Chrien et al. 55(| at Brookhaven 


National Laboratory. Targets of lithium, carbon, oxygen, and aluminum were placed in a 
TT’*' beam at 800 MeV/c and the outgoing proton spectrum was measured. The underlying 


idea behind the experiment can be stated as follows 


56| . If the binding of 77 by a nucleus 
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takes place in the reaction 


7r+ +^A ^p + T]+ ^(A-1) ^ p +^(A-1) , 


then a resonance-like peak will appear in the outgoing proton spectrum. However, the 
expected peak was not observed. Post-analysis has shown that the negative result was due 
mainly to two reasons. First, there was a huge background of proton events. Second, the p 


was produced at high momenta, unfavorab. 
To reduce the background events, Lieb 


o its capture into the Is nuclear orbit. 


57l | has proposed to study the reaction 


TT^ + O p + p + O) ^ p + {tt +p)+ X 

in an experimental setup that favors producing p at rest so as to maximize its capture by 
the nucleus and then take advantage of the A^*(1535) dominance —)■ N* 

7i~ + p) to detect a nearly back-to-back 7r“p pair in coincidence with the outgoing fast 
proton. A peak was observed, but it was located at the area where the detector efficiency is 
limited. Hence, as pointed out in Ref. jd], it is unclear whether the peak was due to the poor 
detector efficiency for which the data were not corrected. Nevertheless, the triple-coincidence 
approach proposed above has since been employed in many experimental studies and has 
successfully reduced background events. 

The search for an p bound state in pion double-charge-exchange (DCX) reaction leading 
to the double isobaric analog state (DIAS), ^®0(7r’'', 7r“)^®Ne(DIAS), was carried out at 


LAMPF by Johnson et al. 5^. At pion energies above the p production threshold, the 


DIAS can be reached via the tt7 r° ^ vr” path as well as via the 7r+ —)■ ^ 7r“ path, 

as illustrated in FigJ2J In FigJS] we show in detail the reaction diagrams of the p-mesic 
nucleus formation. This LAMPF experiment was based on the theoretical calculation 
of the DCX reaction ^^C(7r’'', 7r“)^^0(DIAS). The calculations indicate that when there is 
no formation of bound state of p in the interference between the amplitudes (a) and 
(b) shown in Figj2] will not lead to any new structure in the excitation function of the 
DCX reaction. On the other hand, if there is formation of then the interference 

among the three amplitudes in Fig|2]will produce a resonance-like energy dependence in the 
DCX excitation function near the p threshold. Furthermore, this resonance-like structure 
depends on the pion momentum transfer: greater is the momentum transfer, stronger will be 
the signature. It is reasonable to expect that the excitation functions of the DCX reactions 
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FIG. 2: Reaction diagrams of (a) the vr'*' —^ —>• 7r“ amplitude; (b) the vr'*' —>• r/ —?> 7r“ amplitude; 
(c) the 7r+ —)• T/ —>• 7r“ amplitudes due to bound r]. The shaded, open, and solid circles denote, 
respectively, the GAf^TriV and matrices, and the nuclear vertices. The shaded multiple 

lines denote the 77 -mesic nucleus . 



FIG. 3: Formation process of an 77 -mesic compound nucleus. The decay process [right half of Fig. 
2 (c)] has a similar structure. 
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>'^®Ne(DIAS) and ?-^‘^0(DIAS) have a similar energy dependence. Data analysis 
of the LAMPF experiment showed a visible structure in the predicted energy region of the 


DCX excitation function (Fig.3 of Ref. (58|). However, the statistics of the data was poor. 
Thus, the observed structure was not statistically significant. Experiments having better 
statistics will be greatly valuable. 

A large number of experiments designed to search for ? 7 -mesic nuclei make use of transfer 


reactions. An example was the COSY-GEM collaboration experiment 
the r^-^^Mg bound state in the reaction 

p +^\\ (— )-^He -I- p -|-^^Mg) —)■ ^He -|- vr” -|- p -|- X . 


60[ | designed to detect 


(4.1.1) 


In order to maximize the probability of having the produced p being captured in the Is 
nuclear orbit, the emerging ^He was detected in the forward direction (he., zero degree). 
In this forward geometry the beam momentum is entirely transfered to ^He, leaving the 
produced p at rest. If the p is bound, it cannot emerge as a free p. Instead it interacts with 
a target nucleon and emerges as a pion. For example. 


p + n ^ N*(1535) ^ 71 +p. 


(4.1.2) 


Because the initial p has zero momentum, the emerging 7r“ and p would be back-to-back 
if the neutron had no Fermi motion. With the Fermi-motion, the 7r“ and p will lie in two 
opposite but back-to-back cones. The ^He-7r“-p triple coincidence techniques were employed 
to reduce background events. 

The data are shown in Figj4] where E (denoted BE in Ref. 60(]) represents the real part 
of the p binding energy for events having E < 0. Events having E > 0 correspond to 
an unbound p. The experimental spectrum exhibits a peak structure centered at E = 
— 13.13 ±1.64 MeV with a half-width F/2 ~ 5.1 ± 1.5 MeV. The signihcance of the peak 60| 
is 5.3(j, indicating the existence of an p—‘^^Mg bound state, the mesic nucleus ^^Mg^. 


Q 


Using the factorization approach 1 section llll BD with A = 30 MeV and the BL model 
for pN, we obtained the binding energy Em = —6.5 MeV and half-width F/2 = 7.1 MeV for 
the p—'^^Mg bound state MJ- The large difference between the experimental and theoretical 
values of the binding energy has motivated us to reexamine how the application of a theory 
can take into account the actual set-up of an experiment. 

The usual theoretical approach is to calculate only the following multi-step reaction pro¬ 
cess which we denote as Process M (M for mesic-nucleus formation): 
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FIG. 4: Spectra obtained with (a) potential V giving Em — iT/2 = —(6.5 + i7.1) MeV (dashed 
curve), and (b) potential V giving Em — iT/2 = —(8.0 + 9.8i) MeV (solid curve). The data are 


from Ref. 


60|. 


p + ^^Al rj + +^He 
25Mg^ 

^ 25Mg (tt” + p)+X . 

However, the off-shell rj produced in the intermediate state can also be scattered by the 
residual nucleus and emerge as a pion, without being captured by the nucleus. We denote 
this multi-step reaction process as Process S (S for scattering): 
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p + ^'^Al —)■ ?7 + +^He 

; 

^ 25Mg (tt” + p)+X . 

These two processes are illustrated in FigjSl We emphasize that because these two reaction 
paths lead to the same measured dual state, they cannot be distinguished by the experi¬ 
ment. Consequently, in theoretical analysis one must take coherent summation of the two 
amplitudes to account for the quantum interference between them. We, therefore, ht the 
experimental spectrum by using the sum of two amplitudes: 


® I f S f M I ® 


< k'\V{w)\k>^ 


< P\V{w)\^p >< ^\V{w)\k> 


E — {Efjd — iT/2) 


(4.1.3) 


where V is given by Eq. fl3.2.ip and w = \/l + E. The 'ip is the wave function of bound rj, and 
T is its adjoint (p.l20 of Ref. 30j). We have noted that in the threshold and subthreshold 
regions, //-nucleus interaction is isotropic and that the matrix elements <k'\V\k> are nearly 
constant for k and k' between 0 and 100 MeV/c. Because of these aspects of the //-nucleus 
interaction and the experimental selection of events corresponding to // being produced 
nearly at rest, Eq. (l4.1.3p can be evaluated at \k\ = \k'\ ~ 0. 


In Eq. fl4.1.3p there is only one parameter a and its sole role is to adjust the overall 
magnitude. We emphasize that we used the same micrscopic theory-based V in calculating 
fg and and that the values of Em and r/2 were kept hxed, he., they were not htting 
parameters. Furthermore, we square the sum of the amplitudes, in marked contrast to using 
the sum of a squared background amplitude and a squared Gaussian ampltude. Hence, 
interference effects between the amplitudes are present in our analysis while they were absent 


in the COSY-GEM ht 




Upon introducing the above-mentioned factorization result {Em = —6.5 MeV and 
F/2 = 7.1 MeV) into Eq. (l4.1.3p . we obtained from the ht the overall scale factor a = 4.2 
(counts/fm^) and a spectral distribution peaked at —10.5 MeV. The distribution is shown 
as curve (a) in FigJU The 4.0 MeV downward shift from —6.5 MeV indicates clearly the 
importance of the interference ehect. Ehects of the nuclear medium modihcation of the pN 
interaction, in particular the true pion absorption, was also examined by Haider and Liu 


using the model of Chiang et al. 


62| . They found that the inclusion of true pion absorption 
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(b) 


FIG. 5: (a) Reaction diagram of /g. (b) Reaction diagram of f^. The wavy and multiple lines 
represent, respectively, the r] and ^®Mg. The open oval denotes the ry-nucleus interaction. The 
filled line in (b) denotes the mesic nucleus. 


gave an = —8.0 MeV and T/2 = 9.8 MeV for the M process. Upon introducing these 
latter quantities into Eq. fl4.1.3p . they obtained again the same overall scale factor, a = 4.2 
(counts/fm^), and a spectral distribution peaked at -12.5 MeV (shown as curve (b) in FigJll). 

We have thus seen the importance of interference effect arising from two reaction am¬ 
plitudes. It is worth emphasizing that effects of quantum interference are often crucial 
in understanding the data. For example, researchers were once puzzled by the observed 
“abnormal” cross-section ratios cr[^^C(7r+, 7rN)]/(T[^^C(7r“, ttN)] in the A(1232) region. The 
observed ratios were later well explained when the interference between quasifree and non- 
quasifree reaction amplitudes was taken into account by Ohkubo and Liu 63[ |. 

A different approach to the analysis of the COSY-GEM data on ^^Al was given by Fried¬ 


man et al. 


64l | in which only 


he M Process was considered. By using one of the strongest 


rjN scattering-length models [25| at an appropriate subthreshold energy, they were able to 
obtain for ^®Mg^ binding energies E^d ranging from —14.8 to —19.4 MeV and half-widths 
r/2 between 1.9 and 2.9 MeV. One thus sees that the calculated peaks overestimate the ob- 
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served peak (—13.1 MeV) while the calculated half-widths underestimate the observed one 


by a factor of two. It was argued in Ref. 


64l | that the very small calculated width was due 


to the neglect of true pion-absorption contributions to in-medium rjN amplitude. However, 
no quantitative estimate was given. 

From the view point of nuclear theory an interesting question arises, namely, which 
reaction mechanism is correct? Is it the two-amplitude mechanism leading to cross sections 
proportional to \S + or, is it the one-amplitude M mechanism giving cross sections 
proportional to |Mp ? We believe that the answer lies in obtaining high-statistics data. 
This is because when only the M process is considered, theoretical calculation will give rise 
to a symmetric peak centered at Em with a width F. On the other hand, when both the M 
and S processes are included, the interference between them will result in an asymmetric 
peak (see both curves in FigH]). Current data neither contradict an asymmetric peak nor 
rule out a symmetric one. Future high-statistics data should yield a clear answer. 


65[ | were carried 


Various proton- and deuteron-induced transfer-reaction experiments 
out to search for bound 77 in ^He, ^He, and To date, no bound state has been observed 
in these light nuclei. The photon-induced experiment 


7 He —)■ 7 r° -I- p -I- X 


(4,1,4) 


was performed at MAMI in Mainz 66 (]. A peak structure was seen near the rj threshold 
and was interpreted as an evidence of the mesic nucleus ^He^. However, the analysis of a 
later experiment having much higher statistics revealed that the peak was the result of a 
very complicated structure of the background and, hence, could not support the previous 


conclusions of mesic-nucleus formation 




In view of our discussion on the minimum nuclear mass number needed for the formation 
of an p-mesic nucleus (section [III Bll . we believe that more experiments searching for bound 
states of Tj in nuclei having a mass number A > 12 should be the logical next step. 

B. Final-state-interaction method 

The differential cross section of a two-body to two-body reaction is given by 


^) = (^) I (V’ I V 17+') 


(+)\ |2 


(4.2.1) 
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For nuclear r] production, v is the r^-production potential, p is the beam-target relative 
momentum in the initial channel and k is the ry-nucleus relative momentum in the final 
channel, with and being the corresponding scattering wavefunctions. 


Watson 


54l | showed that in the threshold region, if 


(a) the two-particle scattering in the final channel is dominated by s-wave, 

(b) the primary production (denoted C) is nearly independent of k (apart from energy 
conservation), and 

(c) the interaction between the two final particles is confined in a small region r < a^, such 
that ktts <C 1, 

then the following final-state interaction (FSI) approximation holds: 

I <'I <+’)« i F= ^ 


k (cot S - i)' 

where 6 denotes the s-wave phase shift of the two-particle scattering in the final state. In 
the literature, is termed the enhancement factor. The FSI approach consists in fitting the 
data with the following equation: 




(4.2.3) 


The low-energy expansion of k cot 5 is 


A:cot (5 = - -I- -re/c^ -|- .... (4.2.4) 

a 2 

with a and Vg denoting, respectively, the ? 7 -nucleus s-wave scattering length and effective 
range. In the scattering-length appproximation (SLA), one uses kcot 6 = 1/a. Hence, 
there are three fitting parameters, namely, |C|, Re[a], and Im[a]. In the effective-range 
approximation (ERA), one uses the first two terms of Eq. (14.2.41) . Consequently, there are 
two more fitting parameters, Re[re] and Im[re]. In Appendix [Al we give the inequalities that 
a and/or r^, must satisfy when there is a quasibound state (Eqs. flA.l.lQp and flA.2.14p i. 

Results obtained from fitting the d + p ^ r^-l-^He reaction by different groups are summa¬ 
rized in Table I VI I (For the sake of concise notation, in this subsection the ? 7 -^He scattering 
length and effective range will henceforth be denoted as a and Ve, respectively.) One notes 
from Table IVII that the results of Fits 1 and 2 cannot be used to convincingly determine 
whether there is an r^-^He bound state. This is because the sign of the real part of the 
scattering length are undetermined from FSI fits. Indeed, in the SLA, 











TABLE VI: The s-wave ? 7 -^He scattering lengths, a, and effective ranges, r^, given by different 
fits to various pd —)• r^^He data with the use of Watson’s ESI approximation. The k denotes the 
momenta of p used in the fit. 


Fits 

a [fm] 

re [fm] 

k [MeV/c] 

Refs. 

1 

±(3.8±0.6) -Hi(1.6± 1.1) 


< 75 

m 

2 

±(2.9±2.7) -Hi(3.2± 1.8) 


< 75 

mi 

3 

-2.31 + i2.57 


< 70 

m 

4 

±(10.7 ± 0.8) -Fz(1.5± 2.6) 

(1.9±0.1) ±i(2.1±0.2) 

< 100 


5 

1 4.3 ±0.3 1 ±i(0.5±0.5) 


< 70 

m 


Because in this last equation Re [a] appears as a squared quantity, its sign cannot be deter¬ 
mined. This sign ambiguity also exists in FSI fits using the ERA. 

By defining a = x + iy and re/2 = c + id, we have 


= 


1 -h laTek"^ — ika P /3o + fdik + (32k‘^ + (3sk^ + P^k"^ 


(4.2.6) 


where Pq = l/|ap. Pi = 2yPii, P2 = 1 + {2xc -|- 2yd)PQ, P3 = —2d, and P4 = |rep/4 are 
real-valued coefficients. From the fitted values of Po, Pi, and P3 one can determine |ap, y, 
and d, which, when combined with the P2 given by the fit, will allow one to determine the 
magnitude and the sign of the product 2xc = Re[a]Re[re]. However, one cannot determine 


the sign of Re [a] seperately from the sign of Re [re]. The authors of Ref. 
clear seperation of a and is only possible in a model-dependent way. 


believe that a 


In Fit 3, Wilkin 


19| circumvented the sign ambiguity embedded in Eq. fl4.2.5p with the 


aid of an optical model. First, the sign of Re[a^Ar] was chosen to be positive. The sign-fixed 
QrjN was then used to construct a first-order r^-^He optical potential in the on-shell static 
approximation which, in turn, was used to generate arjSHe = o- The scattering length a 
generated in such a manner has the sign of its real part well-determined by the optical model. 
Finally, the 'K~p —?■ pn and dp —?■ ? 7 ^He data were htted simultaneously by treating both the 
sign-fixed and a as parameters. It was determined that = (0.55 ± 0.20) -|- f0.30 fm 
and a = (—2.31 -|- i2.57) fm. By using the same procedure, it was determined from the 
dd —)■ data that = (—2.21 -f-fl.l). One notes that the r^-^He scattering length, a, 
does not satisfy the bound-state criteria, Eq. flA.l.lOD . in Appendix A-1. Hence, the result 
of Fit 3 indicates that there is no bound state in ^He. On the other hand, the value of 
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given by Fit 3 suggests that an bound state is possible. 


For Fit 4, Mersmann et al. 


53l | only published the value of the scattering length obtained 


in their effective-range approximation to FSI; the effective range in Table I VI I was taken 
from Ref. We have tested the implication of Fit 4 on the existence of r^-^He bound 
state by taking into account all possible combinations of the error bars given in the table. 
We found that only a small number of the combinations satished the existence critera of 
Eq. flA.2.14|) . However, as a whole, the result of Fit 4 could be unreliable because the high rj 
momenta part of the data used in the £t showed substantial p-wave contribution which 
is incompatible with the criteria for using Watson’s FSI theory. 


The scattering length obtained in Fit 5 showed also that the sign o: 
determined by the FSI method alone|68j. It was pointed out in Ref. 


Re [a] could not be 
that there were 


inconsistencies among the experimental data 


43, 1M 


53l |. It was further suggested that one 


should carry out FSI fits by only using data corresponding to rj momentum k < 70 MeV/c 
so as to ensure the s-wave dominance (prerequisite-(a) for using the Watson FSI method). 

In addition to the required s—wave dominance, we recall that prerequisite-(c) for using 
Watson’s method is that the rj momentum, fc, must satisfy kag 1. The maximal 
momentum, kmax, satisfying this inequality can be determined by using the criterion [b 
sm{kmaxCis)/{kmaxtts) — 1- For example, the root-mean-square radius of ^He is ~ 1.22 fm [37|. 
If one assumes 0^=1.2 fm, then sm{kmaxO-s)/(kmaxO's) > 0.97 requires kmax < 0.35 fm“^ ~ 
70 MeV/c. The difference (1 — [sm^kmaxds)/(kmaxCis)]) gives a quantitative estimate of the 
error of the method. 

A general comment regarding the off-shell effects is in order. We see from Eq. fl4.2.2l) that 
the quantity 7F is the on-shell s-wave ? 7 -nucleus scattering amplitude. In other words, the 


FSI method excludes off-shel 
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information. The loss of off-shell information occured when a 


69(] were applied to the original scattering wave function 0*^ ^. 


The on-shell feature of 7F has also been pointed out in Refs. [3|, l68| . We believe that the 


three prerequisites leading to Eq. (I4.2.2j) have minimized the loss of off-shell information. In 
this respect, one could regard the lack of the off-shell effect as a “systematic uncertainty” 
intrinsic to the FSI method. 
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V. SUMMARY AND FUTURE DIRECTION 


The existence of ? 7 -mesic nucleus is a consequence of the attractive interaction between 
7] meson and nucleon in the threshold energy region of the r]N channel. This attactive force 
arises from the strong coupling of the S*!! baryon resonance iV*(1535) to the rjN system. 
While the strength of the attraction is not enough to cause an 7] to be bound on a single 
nucleon, it can cause the r] to be bound on a nucleus, forming a mesic nucleus with a finite 
half-life. The minimum nuclear mass number for forming an r^-mesic nucleus depends on the 
predicted 7]N amplitude. The lightest nucleus onto which an 77 can be bound clearly depends 
sensitively on the real part of the amplitude which is strongly model-dependent. Finding 
the lightest r^-mesic nucleus can, therefore, help differentiate various theoretical models of 
r]N interaction. 

We have shown in Section im that at the formation of an r^-mesic nucleus the t] meson 
interacts with the target nucleon at an energy that is below the 7]N threshold. Consequently, 
only those 7]N models that can be extended to subthreshold energies are relevant to nuclear 
studies of the 77 meson. 

In experimental search for 77 -mesic nucleus, transfer reactions have been frequently em¬ 
ployed. One such reaction has led to the observation of the 77 -mesic nucleus ^^Mg with a 


5.3cr statistical signihcance 


60| . However, searching quasibound 77 -nucleus states in lighter 


nuclei such as ^He, ^He, and has not yet yielded positive results. Searching ? 7 -mesic 
nuclei in medium-mass nuclear systems other than ^^Mg is highly valuable. 

A resonance structure was seen in the excitation function of the pion DCX reaction 
^^ 0 ( 7 r+, 7 r“)^®Ne(DIAS) at energies corresponding to 77 threshold. However, the structure is 
statistically insignihcant. It is, therefore, valuable to repeat this experiment with a much 
higher statistics. If successful, measurement of DCX excitation functions in the same energy 
region for other nuclei, such as "^^Ca and also merit considerations. 

Besides the spectral method, the FSI method has also been used to search for 77 -mesic 
nucleus. However, the FSI method cannot determine the signs of the 77 -nucIeus scattering 
length and effective range without using a theoretical model. Hence, the conclusion of the 
FSI analysis is model dependent. In addition, one must bear in mind that Watson’s FSI 
theory is an approximative theory. The approximative feature is well laid out by the 
three prerequisites of the theory, as outlined in section IV-A. In particular, the dominance 



of final-state s-wave scattering must be ascertained. In other words, the data used in FSI 
analyses must be isotropic in their angular distributions. This isotropy was not fulhlled in 
Fit 4 of Table IVIl We have closely examined all the data used in the hts listed in Table I VI I 
and found that the angular distributions are isotropic only on average (with a large disperson 
of ±5%). Improved data are clearly helpful in future studies. To date, the FSI method has 
been mainly applied to^He. But no convincing evidence of bound state has been found. It 
is equally possible that rj cannot be bound onto a nucleus as light as ^He, as indicated by 
many rjN models. We, therefore, believe that experiments searching for medium and heavy 
rj-mesic nuclei should be given priority in the next step of research. 

Independent of whether an r] can be bound onto a light nucleus, studying rj production off 
a light nuclear system is important. Firstly, because in the production of a physical 77 , the 
basic rjN interaction takes place at energies above the rjN threshold, analysis of rj production 
can, therefore, test an rjN model in an energy domain very different from that relevant to 
the ry-mesic nucleus formation. In addition, the presence of only a few target nucleons in 
a light nuclear system makes calculations of multistep processes involving each individual 
nucleon feasible. Detailed multistep calculations have been reported for pion-induced 0 


and proton-induced 


HEi 


productions. The production data in Ref. 


70| are well described 


by the Bhalerao-Liu model of rjN interaction. Measurements of various differential cross 
sections of r] production and more microscopic analyses of the data are called for. 

The existence of nuclear bound states of the rj in large nuclei creates the opportunity for 
using these mesic nuclei as a laboratory for studying the behavior of an rj meson in a dense 
nuclear environment. This is because the inner region of a medium- or heavy-mass nucleus 
are more close to a dense nuclear medium than the few-nucleon systems are. We mention, for 
example, the suggestion that 77 bound states in nuclei are sensitive to the singlet component 
in rj and can be used as a probe of flavor-singlet dynamics [^. Hence, 77 -mesic nuclei can 
improve our understanding on the 77 - 77 ' mixing. There is also theoretical work indicating 
that a dense nuclear medium can have very different effects on the hyperon A(1405) and the 
baryon V*(1535) j^. These predicted effects can be checked by means of 77 -mesic nuclei. 
There are also suggestions that the formation spectra of the 77 -mesic nuclei can be used to 
check the prediction by chiral doublet theory on the mass difference betweem the nucleon 


and the iV*(1535) in a nuclear medium |75|,176(. 


We would like to conclude this review by emphasizing the importance of searching 77 - 
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nucleus bound states in medium and heavy mass nuclei. We believe that the successful 
observation of the encouraging structure in the excitation function of pion DCX 

reaction, the progress made in FSI studies, and the mastering of triple-coincidence and 
7] —)■ 77 detection techniques have laid down a solid foundation for future successful searches 
of ? 7 -mesic nuclei. 
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Appendix A: Analytical Relations between bound-state and scattering observables 

The s-wave scattering amplitude is given by 

where S{k) = (cot <5 -|- i)/(cot S — i) is the S'—matrix, S is the phase shift, and k is the c.m. 
momentum. For potentials that are exponentially bound, one has the following low-energy 
expansion; 

kcot6 =— I — r^k'^ + .... (A2) 

a 2 

where a denotes the s-wave scattering length and the s-wave effective range. When an 
optical potential is used in the calculation, the quantities 6, a, and are all complex-valued. 
If there is a bound state (also termed quasibound state) of complex momentum kpoi , then 
the S'-matrix has a pole at kpoi. It follows that in Eq. flAlj) 

kpoi cot 6 -ikpoi = 0 . (A3) 

In this appendix, we derive the interaction-model independent analytical relations between 
the binding energy, width, kpoi, scattering length a, and the effective range r^. 
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1. The scattering length approximation 


Equation flA2p shows that the hrst term dominates when k is very small. Hence, one may 
approximate the low-energy expansion by using 


k cot 8 = — . 
a 


Equation flASp then gives kpoi = —i/a and 


(A.1.1) 


p ^_1 

^pol 2 * 

In the complex a-plane, we may express the scattering length by 


(A.1.2) 


a = |a| exp(i 7 ) = x + iy , 


(A.1.3) 


with 


Hence, 


X = Re[a] = |a| cosy, y = Im[a] = |a| siny, y = arctan f—j . (A.1.4) 

\ X / 


p = _^ 

^pol |„|2 


exp(—2iy) . 


The complex energy B is, therefore, given by 

H = ^ = -\B\ exp(-2^y) , 

2/i 

where y is the reduced mass of the bound particle and 

1 


B = 


2y\a\ 


In the Cartesian representation. 


r 


B = E — i— = u + iv , 


(A.1.5) 


(A.1.6) 


(A,l,7) 


(A.1.8) 


where E {E < 0) and r/2 (T > 0) denote, respectively, the binding energy and half-width 
of the bound state. It follows from Eqs. flA.l.Sp . (lA.l.bp . and flA.l.Sp that 


u = E = 


— y'^ 

2p|a|^ ’ 


(A.1.9) 


_ r 2xy 


2 2y\a\'^ 
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(A.1.10) 














Because u < 0 and u < 0, in what follows we will often write u 
whenever it is more convenient. 

If we denote 

^pol R 'll 5 

then 

kl^i = R^ + 2iRI , 


u\ and V = —\v 


(A.l.ll) 


(A.1.12) 


and, from Eqs. flA.l.hp and flA.l.Sp . 

2/iE= R^-P = 


cos 27 , 


(A.1.13) 


2/i ( - ) = -2RI = -7777 sin 27 . 


(A.1.14) 


Because E < 0 and B > 0, Eqs. flA.l.lSp and flA.l.ldp are, respectively, equivalent to 


{R^ — P) <Q and RI <Q . 


(A.1.15) 


The hrst inequality requires |i?| < |/|. In addition, a decaying outgoing wave of a bound 
state requires / > 0. The second inequality then leads to i? < 0. In summary, a bound-state 
requires simultaneously 

|i?| < |/|, / > 0, i? < 0 , (A.1.16) 

which indicate that the bound-state poles kpoi are situated in the second quadrant (but 
above the diagonal line) of the complex fcpoi^plane. 

Indeed, upon solving Eqs. (lA.1.13p and (IA.1.1411 for R and /, one has 

R = — \J-\- pV? -{- = —\J /U.(—|m| -|- pP -f vP , (A.1.17) 


I = 


\Jn{—u + pp + p \u\ + pp +P^) . 


(A.1.18) 


In choosing the branch of the square roots, the properties associated with the physical 
domains discussed above have been used. The R and I clearly satisfy all the three conditions 
stated in Eq. flA.1.161) . 

Equations flA.1.13p and flA.1.14p further indicate that E < 0 and T > 0 require, respec¬ 
tively, that cos 27 > 0 and sin 27 < 0. This in turn requires 37 r /4 < 7 < vr. The complex 
scattering length, a, is therefore situated in the second quadrant but below the diagonal 
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FIG. 6: The complex scattering length plane. The physical domain is the entire lower triangular 
region of the 2”'^ quadrant. The meaning of the solid circles and the crosses are given in the text. 


line in the complex a-plane. In other words, in the a—plane the scattering length satishes 
simultaneously 


Im[a] > 0, Re[a] < 0, |Im[a]| < |Re[a]| . 


(A.1.19) 


The third inequality was hrst given in Ref. Mi- The physical domains in the a-, kpoi-^ and 
R-planes are shown in FigslGHH], respectively. When the polar angle, 7 , in the a-plane turns 
counter clockwise, the corresponding polar angles in the kpoi- and R-planes turns in the 
opposite direction. 


We now proceed to express the scattering length a in terms of the binding energy B. By 
using of Eq. flA.1.7D . we can rewrite Eqs. flA.l.Qj) and flA.l.lOj) as 


u =-2ij.\B\\x^ - y^) , 


(A. 1 . 20 ) 


V = 2y\B\‘^{2xy) , 


(A. 1 . 21 ) 


where, because of Eq. flA.1.8|) . \B\ = \/v? + v^. The inverse mapping, R —)■ a, is obtained by 
solving the above coupled equations and the result is: 


X = 


( 2\B\^ 



+ ^/v? + v'^ 


I 

mi^ ’ 


(A. 1 . 22 ) 
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FIG. 7: The complex kpd plane. The physical domain is the entire upper trangular region of the 
2nd quadrant. The arc near the origin indicates the corresponding polar-angle range. 



FIG. 8: The complex energy plane. The physical domain is the entire 3'’'^ quadrant. The arc 
indicates the polar-angle range. The trajectories shown by the downward dashed line and the 
horizontal linked-dotted line are explained in the text 
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y = 


(2^) \/-M + '/^^ = 


R 


(A.1.23) 


2\B\^J V i-i ■ • 

In the inverse mapping the object trajectory is on the i?-plane and the image trajectory 
is on the a-plane. We have considered the following two cases as examples. Case I corre¬ 
sponds to hxing E but varying r/2 in the direction given by the downward arrow along the 
dashed line in FigJHl Case II corresponds to hxing the value of r/2 while varying E along 
the horizontal linked-dotted line in the direction shown by the leftward arrow in FigjHl The 
resulting a (calculated with 2p = 5 fm“^) are shown, respectively, as the left-to-right solid 
circles and the descending crosses in FiglHl One notes that the polar angles, 7 , of the suc¬ 
cessive left-to-right solid circles in Fig J 6 ] turn clockwise while the polar angles of the original 
5-points along the downward dashed line in Fig I 8 ] turn in a counterclockwise direction. This 
reversal of the sense of the turnung is a consequence of the opposite signs in front of the 
polar angles in Eqs. (IA.1.3j) and flA.l.bl) . The above inverse mappling explains the general 
feature given by detailed calculations in Ref. [77 1. 

Discussion on the role of nuclear mass in binding a particle'. From Eqs. flA.1.22p and 


flA.1.23p . one notes that the magnitudes of the real and imaginary parts of the scatter¬ 
ing length a are inversely proportional to ^JJi. For 7 -mesic nucleus, /i is the reduced mass 
of ? 7 , which increases as nuclear mass increases. This in turn indicates that for rj to have a 
given binding energy \E\ in a lighter nucleus, it requires a larger |Re[a]| than that for rj to 
have the same binding energy |5| in heavier nuclei. 

Finally, it is worth pointing out that in the literature, another sign convention of the 
scattering length is also used, namely, lim^^o kcotS = —1/a. We showed in Ref.[78[ that 
this sign difference does not alter the obtained results. 


2. The effective range approximation 


In the effective-range approximation to the low-energy expansion, both the hrst and 
second terms of Eq. (IA2p are retained and Eq. flA3p becomes 

1 


T ikj)oi 0 


(A.2.1) 
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There are three variables, a, rg, and kpoi in this equation. Conseqently, one can only express 
one of the three variables as a function of the other two. By writing 


— = c + id , 

one obtains readily from Eq. flA.2.ip 

{\u\c — \v\d) — I / (2/i) 

— {\v\c \u\d) — R/ (2/i) 

F(c, d) 

where 


(A.2.2) 


(A.2.3) 

(A.2.4) 


F(c, d) = (2/i)|F|"(F + F) + \B\ + 2c{R\v\ - I\u\) + 2d{R\u\ + I\v\) . 


(A.2.5) 


It is easy to verify that when c = d = 0, Eq. flA.2.5p becomes F(0, 0) = |F|, and Eqs. flA.2.3p 
and flA.2.4D reduce to Eqs. flA.1.221) and flA.1.23D . respectively. 

By treating as the unknown in Eq. flA.2.1D . we have the relations 

v{R + ///|ap) — u{I + x/\a\‘^) 


c = 


2id{u^ + 


(A.2.6) 


d = 


v{I + x/lap) + u{R + y/\a\^) 


(A.2.7) 


2iJi{v? + 

Finally, we solve for kpoi when a and are known. As one can see, Eq. flA.2.ip yields two 
solutions for kpoi- However, only the following one is physically meaningful, namely. 



(A.2.8) 

This is because kpoi —i/a as 0, so that the scattering-length approximation is 

recovered. It follows from Eq. (IA.2.8p that 



(A.2.9) 


If |2re/a| < 1, we can make a Taylor’s expansion of the square root. To the order of (2re/a)^, 
Eqs. flA. 2. 8p and flA.2.9p are, respectively, equal to 


kpoi 




(A.2.10) 
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p = 

^pol 


r^ - a 


The condition Re[kp^i] < 0 is then given by 


Re 


fre-a\ 

1 

1 ) 

|a|6 


= -j—^Re[a*^(re — a)] < 0 , 


(A.2.11) 


(A.2.12) 


or 


Re[a^(r* - a*)] < 0 . (A.2.13) 

This last inequality, resulting from the use of limited Taylor expansion, was hrst given in 


Ref. 


68l |. Consequently, in the effective-range approximation the necessary conditions for 


having a bound state are: 


Re[a'^(rg — a*)] <0; / > 0; R < 0 , 


(A.2.14) 


where I and R are obtained from taking, respectively, the imaginary and real parts of 
Eq. flA.2.10p . When |2re/a| is very large, either the Taylor’s expansion cannot be done or it 
brings no advantage. In this latter case, we can directly calculate Eq. flA.2.8p by using the 
polar representation. 

By successively dehning 


T 

P = i = W\ ? = —= |J|e‘>‘, C=l + { = IC|e‘* 

^ CL 


we can rewrite Eq. (]A.2.8p as 

kpoi 

which gives 


2|p| 


^—iX 




R = Re[kpoi\ = ^ sin A + |C| ^ sin(^ - A)^ , 

I = lm[kpoi\ = ^ cos A - ICI^ cos(^ - A) 

From Eqs. flA.2.16p - flA.2.18p . one readily obtains 


R^-I^ = 


4|p|^ 


-cos(2A) — \(\cos{ip — 2A) -I- 2\(^\'^cos(^ — 2A) 


(A.2.15) 

(A.2.16) 

(A.2.17) 

(A.2.18) 

(A.2.19) 


Again, for a bound state to exist, the quantities R^ — R, R, and / must satisfy the three 
inequalities of Eq. flA.l.lbp . 
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In summary, all the analytic expressions derived in this appendix are interaction-model 
independent as long as the potential of the particle-target interaction belongs to the class 
of exponentially bound potentials so that the low-energy expansion, Eq. flA2p . can be made. 
The only kinematic approximation used in our derivation is kp^i/2fi ~ ^— ji which 
is a very good approximation for bound state problems. We emphasize that the model 
depedence of the interaction dynamics does come into play when one theoretically calculates 
the scattering length, a, the effective range, re, and the binding energies, B. In the effective- 
range approximation when two of these three quantities are calculated (or measured), the 
remaining one will be fixed by the analytic relations. Similarly, in the scattering-length 
approximation, when a (or B) is calculated (or measured), the other variable B (or a) will 
be hxed by the analytic relation. In this respect, these analytic relations can be used to 
check the consistency of calculations (or measurements). 
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